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Abstract 

Within a nuclear model that the deuteron has NN and ttNN components, we derive convolution 
formula for investigating the Drell-Yan process in proton-deuteron (pd) reactions. The contribution 
from the irNN component is expressed in terms of a pion momentum distribution that depends 
sensitively on the ttNN form factor. With a ttNN form factor determined by fitting the ttN 
scattering data up to invariant mass W = 1.3 GeV, we find that the pion-exchange and nucleon 
Fermi-motion effects can change significantly the ratios between the proton-deuteron and proton- 
proton Drell-Yan cross sections, R p d/ PP = & pd / (2a pp ), in the region where the partons emitted from 
the target deuteron are in the Bjorken X2 > 0.4 region. The calculated ratios R p d/ pp at 800 GeV 
" agree with the available data. Predictions at 120 GeV for analyzing the forthcoming data from 

Fermilab are presented. 



PACS numbers: 13.85. Qk, 13.60.Hb, 14.20.Dh, 13.75. Cs 
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I. INTRODUCTION 



Since the asymmetry between the anti-up (u) and anti-down (d) quark distributions in the 
proton was revealed by the New Muon Collaboration [l[ (NMC), a series of experiments 0- 
[i[ on the di-muons (a + u~) production from the Drell-Yan [gj (DY) processes in pp and pd 
collisions had been performed at Fermi National Accelerator Laboratory (Fermilab). The 
objective was to extract the d/u ratio of the parton distribution functions (PDFs) in the 
proton. The information from these experiments and the measurements [j], 0, S| of deep 
inelastic scattering (DIS) of leptons from the nucleon have confirmed the NMC's finding, 
d/u > 1, only in the region of low Bjorken x < 0.35. 

The ratio d/u > 1 signals the nonperturbative nature of the sea of the proton. Its 
dynamical origins have been investigated [9-23] rather extensively. Precise experimental 
determination of d/u for higher x > 0.35 is needed to distinguish more decisively these 
models and to develop a deeper understanding of the the sea of the proton. This information 



will soon become available from a forthcoming experiment |24| at Fermilab. 

In analyzing the DY data on the deuteron [H-|5| and nuclei [25-28||, it is common to 
neglect the nuclear effects that are known to be important in analyzing the DIS data. It is 
well established that the nuclear effect due to the nucleon Fermi motion (FM) can influence 
significantly the DIS cross sections, in particular in the large x region. It is also known 
that the contributions from the virtual pions in nuclei must be considered for a quantitative 
understanding of the parton distributions in nuclear medium. Thus it is interesting and 
also important to develop an approach to investigate these two nuclear effects on the pd DY 
process. This is the main purpose of this work. We will apply our formula to analyze the 
available data at 800 GeV jBf and make predictions for the forthcoming experiment [24|. 

It is instructive to describe here how the DY data were analyzed, as described in, for 
example, Ref. (Hf. The leading-order DY cross sections from pN collision with N = p 
(proton), n (neutron) is written as 

lh " X l7T '' 2 J2 i2 M X ^fN^) + #(*l)/&(* 2 )] , (1) 



dx x dx 2 9M 2 



where the sum is over all quark flavors, e q is the quark charge, f^{x) is the parton distribution 
of parton q in hadron N, and M is the virtual photon or di-lepton mass. Here X\ and x 2 
are the Bjorken-x of partons from the beam (p) and target (N), respectively (see Sec. IV Al 
for explicit definitions of X\ and x 2 ). The DY cross section for pd is taken to be 

da pd da pp da pn 

V 1 ) 



dx\dx2 dx\dx2 dx\dx 2 

Obviously, Eq. (T5]) does not account for the nucleon Fermi-motion and pion-exchange effects. 
To make progress, it is necessary to investigate under what assumptions Eqs. ([1]) and (T5]) can 
be derived from a formulation within which these two nuclear effects are properly accounted 
for. 

We start with a nuclear model within which the deuteron wave function has NN and ttNN 



components. Such a model was developed in the study of nNN system 29] . We will derive 
convolution formula to express the DY cross section in terms of the momentum distributions 
p(p) calculated from the NN component and p(k v ) from nNN component. Since 7rNN 
component is much weaker, it is a good approximation to use the NN component generated 
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FIG. 1. DY process. 



from the available realistic NN potentials [30] • In the same leading order approximation, 
the resulting p(k 7T ) depend sensitively on the nNN form factor. An essential feature of our 
approach is to determine this form factor from fitting the ttN scattering data. This provides 
an empirical constraint on our predictions of the pion effects on the proton-deuteron DY 
cross sections in the un-explored large x region. 

To see clearly the content of our approach, we will give a rather elementary derivation of 
our formula with all approximations specified explicitly. In Sec. Ull we start with the general 
covariant form of the DY cross section and indicate the procedures needed to obtain the well 
known qq — > p + p~ cross section a qq . The same procedures are then used to derive formula 
for calculating the pp and pn DY cross sections from a qq and a properly defined PDFs 
of the nucleon. 

In Sec. we use the impulse approximation to derive the formula for calculating pd DY 
cross sections from a qq , f^, and the momentum distributions Pn{p) for nucleon and p^k) 
for pions in the deuteron. The calculations of these two momentum distributions within the 
considered ttNN model are explained in Sec. IIV1 

In Sec. [V] we develop the procedures for applying the developed formula to perform nu- 



merical calculations of pp and pd DY cross sections using the available PDFs |3lH35| and 
realistic deuteron wave functions j30| . In Sec. IVTj we present results to compare with the 
available data at 800 GeV and make predictions for analyzing the forthcoming experi- 
ment 24j. A summary is given in Sec. IVHI 



II. FORMULA FOR DY CROSS SECTIONS 



The formula presented in this section are derived from using the Bjorken-Drell [36] con 



ventions for the Dirac matrices and the field operators for leptons, nucleons, pions, and 
photons. We choose the normalization that the plane- wave state \k) is normalized as 
(k\k') = 5(k — k') and the bound states of composite particles, nucleons or nuclei, 
are normalized as (& a \$p) = 5^. To simplify the presentation, spin indices are suppressed; 
i.e. \k a ) represents \k a , \ a ) for a particle a with helicity A a . Thus the formula presented here 
are only for the spin averaged cross sections which are our focus in this paper. 

We consider the di-muons production from the DY processes of hadron (Zi)-hadron (T) 
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collisions: 

h(p h ) + T(p T ) ->• fi + (k + ) + iT(k-) + X h (p Xh ) + X T {px T ), (3) 

where Xh and Xt are the undetected fragments, and the four-momentum of each particle 
is given within the parenthesis. In terms of the partonic qq — > 7 — > fi + ^~ mechanism, 
illustrated in Fig. [TJ the covariant form of the di-muons production cross section can be 
written as 



d*= AT ,_ _ h - a _2ii/»^o^; ^H^ fc -)^Pr.g). ( 4 ) 



(2?r) 4 1 dk+ dk. 1 

4[(p h • p T ) 2 - mlml} 1 ' 2 (2vr) 6 2E+ 2E_ J 4 

where and my are the masses for h and T, respectively, Ey- = \Jk± + are the energies 

of muons fi , and q = k + + k_ is the momentum of the virtual photon. The leptonic tensor 
is defined by 

r(k + ,k_) = (2n)%2E + )(2E_)(kJ^fmO)(0\fmk + k-). (5) 
Here taking summation of lepton spins is implied. The leptonic current is 

f{x) = e^(x) 7 ^(x) , (6) 

where ipn(x) is the field operator for muon, and e = v 'Ana with a = 1/137. By using the 
definitions Eqs. (JSJ) and fl6]), it is straightforward to get the following analytic form of the 
lepton tensor 

r{k + , fc_) = -4e 2 [k$k v _ + k\k» - <r(k + ■ fc_ + ml)]. (7) 

Within the parton model, the hadronic tensor in Eq. (j3J) is determined by the current 
Jfi{x) carried by partons q or q 

F^(Ph,Pr,q) = ^2 i 27l f(2E h )(2E T ) J dp Xh dpx T S 4 (Ph+PT-Px h -Px T -q) 

x (PhPr 1^(0) \p Xh dpx T ) {Px T Px h \J V (0) \PhPt)- (8) 

Here it is noted that, throughout this paper, we shall take the summation (average) for the 
spins of final (initial) particles appearing in hadronic tensors. 

J A X ) = ^^(xh^Or) , (9) 
g 

where ip q (x) is the field operator for a quark q with charge e q e; i.e e u = 2/3 and id = —1/3 
for up and down quarks, respectively. 

The above covariant expressions are convenient for deriving the formula that can express 
the hadron-hadron DY cross sections in terms of the elementary partonic qq — > fi + fi~ cross 
sections. To get such formula, we first show how the elementary qq — > n + fi~ cross section can 
be derived from Eq. (j3J) with h = q and T — q. We then derive the formula for calculating 
proton-nucleon DY cross section. 
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A. qq — > fi + fi cross section 

Explicitly, Eq. (j3J for the q{p q ) + q(p q ) — > n + {k + ) + process is 

The next step is to replace the intermediate states \px h Px T ) by the the vacuum state |0) in 
evaluating the hadronic tensor Eq. ([8]). We thus have 

F$(p 9 ,P 9 ,q) = (27i)\2E q )(2E,)(p,p g \J^)(0\J u mP q ,P^\p q +P, ~ ?) • 

(11) 



Substituting parton current (J9J) into Eq. ffTTj) . the hadronic tensor Fffi then has a form that 
is the same as the leptonic tensor f^ v defined by Eqs. (jSJ) and ([6]) except that the momentum 
variables and charges are different. By appropriately changing the momentum variables in 
Eq. (j7|), we obtain 

F$(j> q , Pq ,q) = ^e) 2 [py,+py,-g^(p q - Pq + m 2 q )}5 4 (p q +p,-q). (12) 

Here, compared with the lepton tensor case [Eq. ([7j)], the difference of factor 4 is because 
the average of quark and anti-quark spins is taken for this case. By using Eqs. ([7]) and ffT2l . 
Eq. ffTOj) for the cross sections of q(p g ) + q{p q ) — > fi + (k + ) + n~{kJ) can then be written as 



7 no/ \ (2tt) 4 1 dk + dk- 1 . 4 . . 

d<J ijPqi Pg) = ^ )2 _ /2 — — — (p, + p, - g) 



x8e 2 e 4 



£;+ • p q k_ -p q + k_- p q k + ■ p q + m h m 



(k + + k_f t _ 2 ( Pq +Pg) 2 

> 2 



.(13) 



It is convenient to express the qq DY cross section in terms of the invariant function 
<? 2 = (Pq + Pq) 2 = (k + + k_) 2 . After some derivations and accounting for the color degrees 
of freedom of quarks, we arrive 



g 2 6g 3N c [g 2 -4m 2 ]V2 ^ g 2 J ^ 



,2 



" ^(g 2 -(p g + p g -) 2 ),(l4) 



dg 2 g 2 q 3N c [q 2 - Am 2 ] 1 ' 2 \ g 2 J V g 



,2 



where iV c is the number of colors. Taking iV c = 3 and considering g 2 ^> m 2 and g 2 3> m 2 
we then obtain the familiar form 



da qq (Pq,Pq) 47ra 2 „ 2 , 2 2 \ 

^ = "9^2- e ^ V ~ (P« + Pff) ) • ( 15 ) 

The above expression is identical to the commonly used expression, as given, for example, 
by the CETEQ group 

In Sec. Ill Bl we will derive formula expressing the pN cross sections in terms of da qq (p q , p q ) / dq 2 
given in Eq. f lT5|) . 
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B. pN DY cross sections 

To simplify the presentation, we only present formula for q in the projectile p and q in the 
target N. The term from the interchange q -H- q will be included only in the final expressions 
for calculations. 

Equation (JU) for the p(p p ) + N(p N ) ->■ p + (k + ) + fi~(kJ) + X p (p Xp ) + X N (p XN ) process is 

= iis^fe^ t^fefe^-^^^'l (16) 

where the hadronic tensor, defined by Eq. (JHJ), is 

F^( Pp ,p N ,q) = (27r) 6 (2E p )(2E N ) £ J d p Xp dp x j\p p + p N - p Xp - p XN - q) 

x P ,x N 

X (PTVPp |^(0) ) {PX N PX P \J„(0) \PpPN ) • (17) 

Within the parton model, the DY cross sections are calculated from the matrix element 
(qq\ ,7^(0)10) (0\J v (0)\qq) which is due to the annihilation of a q (q) from the projectile p and 
a q (q) from the target N into a photon. To identify such matrix elements, we insert a 
complete set of qq states (omitting spin indices) 



dPqdpq\p q Pq){PqP q \ 

into Eq. ( TPTj) . We then have 

F^(p p ,p N ,q) = (27r)%2E p )(2E N ) ]T J dp Xp dp x j\p p + p N - p Xp - p Xpf - q) 

Xp,Xfir 

X J dp q dpq J dp^^p^fx^^NW^X^^X^q^^X^qWN) 

X (PqPq | J M (0) |0> (0| J„(0) IPX) • (18) 

By momentum conservation, the overlap functions in the above equation can be written 

as 

{PqPx\p P ) = (Px\bpr q \p P ) 

= <Pp P (Pq,Px)S(Pp-Pq-px), (19) 

where bp q is the annihilation operator of quark q. 

By using the definition (1T91 . Eq. (ITS]) can be written as 

F^(p p ,p N ,q) = (27r) 6 (2E p )(2E N ) ]T J dp Xp dp x j\p p + p N - p Xp - p XN - q) 

Xp,Xn 
|2 



X y ^|0p p (Pg,PX p )| $(Pp-Pq-PX p ) 
X / dpq\(pp N (Pq,p XN )\ 2 S(p N - Pq - f XlK 



x(PqP q \U0)\0}(0\Ju(0)\P^)- (20) 
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The evaluation of Eq. (|20|) needs rather detail information about the undetected fragments 
X p and Xn because of the dependence of 8 A {p p + pn — Px p — Px N — q) on their energies p Xp 
and p°x N - To simplify the calculation, we follow the common practice to neglect the explicit 
dependence of the energy p x and p x of the undetected fragments. This amounts to using 
the approximation p x ~ e\ and p\ N ~ e 2 , where E\ and e 2 are some constant energies, to 
write 



$ 4 (Pp +Pn~ Px v - Px N ~ q) ~ 5(p p + p N - p Xv - Px N - q)5(p° p + p° N - e 1 - e 2 - q°). 



(21) 



We now define 



= J2J d Px P \(t)p P {PQ^x v )\ 2 5{p p -p q -p Xp ), (22) 

Xp 

for the projectile p, and 

/|v (P?) = Yl J d P x " I fe» ) 1 25 ^ -Pv-Pxn), (23) 

for the target AT. These two definitions and the approximation (12T|) allow us to cast Eq. (j20|) 
into the following form 

F^(p P ,PN,q) = (2n)%2E p )(2E N )J2 [ dp q dp q f^ 
x(^|J M (0)|0)(0|J,(0)|^) 

x5(p 9 + p 5 - g)«J(p; + - e l - e 2 - g°). (24) 

We next make a reasonable approximation that ei (e 2 ) in Eq. (|24|) is the difference between 
the energy of the projectile p (target N) and the removed parton q (q); namely assuming 

5(p° p + p° N -e 1 -e 2 -q°) = 5 ((p° p - Cl ) + (p° N - e 2 ) - g°) 
Then Eq. (I24I) can be written as 



~*(p5+p5-9°)- (25) 



F^( Pp ,p N ,q) = £ / dp q dp q fl(p q )fl N (p 



^ E P E N 

E q Eq 

x{(27r) 6 (2£ 9 )(2i^(2^ 

(26) 

The quantity within the bracket {...} in the above equation is just the hadronic tensor 
F^(p q ,p q ), defined in Eq. ffTTj) . for gg system. We thus have 

*ff(Pp,PN, q) = E/ ^^(^)4(^)§§^^5(^^ ? -) ■ (27) 
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FIG. 2. Impulse approximation of pd DY process. 



Substituting Eq. ( 1271) into Eq. ( fT6l) . we then have 

,E p E N 



1 dk-L. dk- 1 



(2tt)< 



(2tt) 6 2E+ 2E_ q 4 



E q E q A[(p p -p N ) 2 
r(k + ,k^Ffi(p q ,p,) 



in 



(28) 



The quantity in the bracket {...} of the above equation is precisely what is in the bracket 
{...} of Eq. ( [TO]) for the qq — > / u + yT process. Accounting for the difference in flux factors 
and extending Eq. (128]) to include the q -<->• q interchange term, the full expression of the pN 
DY process is 



da pN (p p ,p N ) 
dq 2 



X 



[(p« • p?) 



411/2 



E p E N da^ip g ,pg) 



[(Pp ■ Pn) 2 - mjmjf)} 1 / 2 E q Eq dq 2 



(29) 



where da qq (p q ,Pq)/dq 2 is the qq DY cross section, as defined by Eq. ( [15]) . 

We now examine the physical meaning of the functions fi p (p q ) and f^ N {Pq) in Eq. (129]) . 
By using the definitions (fl9|) for (fip p {pq,Px p ) and Eq. ( 122"]) for fp p (p q ), it is straightforward 
to show that 



(PpI & p- & pJPp) 
(PpIPp) 



(30) 



Thus fp (p q ) is just the probability of finding a quark q with momentum p ? in a nucleon state 
Note that this simple interpretation of fp(p q ) is due to the use of the approximations 
Eqs. (j2~T]) and (125]) . If we depart from these two simplifications, we then need the spectral 
function of the nucleon in terms of parton degrees of freedom to make calculation for DY 
cross sections. Such information is not available at the present time. 
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III. PROTON-DEUTERON DY CROSS SECTION 



In this section, we derive formula to express the proton-deuteron (pd) DY cross sections 
in terms of da qq (p q ,Pg)/dq 2 of Eq. ffl~5l) for the elementary q + q — >• fi + + process. To 
simplify the presentation, we only explain the derivation of the formula for a quark q emitted 
from the projectile p and an anti-quark q from the target d. The terms from the interchange 
q <H- q will be included only in the final expression of the cross section. 

For the p{p p ) + d(pd) — > H + {k + ) + + X p + X d DY precess, Eq. (jU) gives 

^ = 4[( Pp ■ ^L^n ^W-A^' k - ) ^' P " qy (31) 

To proceed, we need to define a model for generating the deuteron wave function. Here we 
follow the nNN studies 29|, |37| to consider a nuclear model within which the deuteron wave 
function has two components 

\*d) = \$d) + I&ttNn), (32) 
where ^ is the NN component. In the following two subsections, we derive formulas for 
calculating the contribution from each component of the deuteron wave function to the 
pd DY cross sections. 



A. Contribution from nucleons 

We assume that the pd DY process takes place on each of the nucleons in the deuteron, 
as illustrated in Fig. |2j In this impulse approximation, the hadronic tensor for a deuteron 
target can be obtained by simply extending Eq. ( TT7|) for the pN to include a spectator 
nucleon state \p s ) in the sum over the final hadronic states. We thus have 

F Jlu(PpiPd, q) = {2itf{2E p )(2E d ) ^ J d Psdp Xp dpx N 5 A (p p + p d - p Xp ~ Px N - p a - q) 

x (^p d PpM°)\Px P Px N Ps)(PsPx N Px p \Ju{0)\^ , (33) 

where \& Pd ) is the NN component of a deuteron moving with a momentum p d . We expand 
\$ Pd ) in terms of iVW plane- wave states 

\®p d ) = J dp N dp2$ Pd (pN)S(Pd -Pn- P2)\pnP2)- (34) 

Keeping only the contributions due to a parton in \pn) of the above expansion of <& Pd and a 
parton from projectile state \p p ), the current matrix element in Eq. fl33|) becomes 



(% d P P \M°)\Px P Px N Ps) = j dp N $* Pd (p N )8(pd- Pn - Ps){pnP p \J»(0)\px p Px n )- (35) 
By using Eq. fl35|) . Eq. fl33l can then be written as 

x{(2n) 6 (2E p )(2E N ) 

x / dpx p dpx N 5 A (p P + p N - Px p ~ Px N ~ q) 

v v__ J 



Xry.X 



x {PnP p 1^(0) \p Xp Px N ) (px N Px P \J V (0) \p p Pn)} ■ (36) 
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We see that the quantity within the bracket {...} in the above equation is identical to 
F£»(Pp,Pn, q) of Eq. (JTTJ) . We then have 

Fff( Pp , Pd , q) = J dp N p Pd (p N )^-F^(p p ,p N , q), (37) 

where 

PpApn) = \%Apn)\ 2 . (38) 



By using Eq. (|34|) . one can show that 

PpApn) = (% d \bl N b fN \% d ), (39) 

where £>j- is the creation operator for a nucleon with momentum p^. Thus p Pc1 {pn) is the 
nucleon momentum distribution in a moving deuteron with momentum p d . We will present 
formula for calculating p P(1 ( P n) in Sec. [TV] 
By using Eq. (13"?]) . Eq. (1HT|) becomes 

^ ^ ( 2?r ) 4 sr \ E v E * 



4[(ft, -^) 2 - m> 2 ]V2 ^ J —^' Ep E N 



N=p,n 

1 c£AL_ dk 



(2w) 6 2E + 2E_ g 4 ' 



\r(k + ,k-)F^(p p ,p N ,q) . (40) 



The quantity within the bracket {...} of the above equation is exactly what is in the bracket 
{...} of Eq. ( fl6l) . Accounting for the difference in flux factor, we obviously can write 

dgP d (p p ,p d ) _ x ^ f _ [(Pp ■ Pn? ~ mlml] 1 ' 2 E p E d da pN (p p , p N ) 

dq 2 ~ J^J dpNp ^ N) [(p p • p d ) 2 — mpWij] 1 / 2 E p En ' (41) 

where da pN {p p ,p^) / dq 2 are given in Eq. ( 12 9 1) . 
Substituting Eq. Q29]) into Eq. (HI]) , we have 

dv pd (p P , Pd ) _ ^ ,., r /,.„ [GVP*) 2 -™J] 1/2 ^ P ^d 



dq 



x If I &)& (ft) + 4r (ft)/| p (ft)] ir ff; • (42) 



B. Contribution from pions 

In the impulse approximation, the contribution from the ttNN component of Eq. (132]) 
to the pd DY cross sections can be derived by using the similar procedures detailed in the 
previous subsection. We find that the resulting formula can be obtained from Eq. (142 j) by 
changing the momentum distribution and parton distributions for the nucleon to those for 
the pion. Explicitly, we have 

da pd (p p ,p d ) _ f, r ,r^[~^ [(ft" ft) 2 -™'] 172 E p E d 



dq 2 



jt ft Kft'ft) ~KJ V E P E, 

x[/| p (ft)/£(ft) + /J (ft)4(ft)]^r^ , (43) 
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where /| and f$ are the PDFs for the pion, and the pion momentum distribution in a 
moving deuteron with momentum p d is defined by 

P Pd {K) = {^NN. Pd W^HM^NN,p d )- (44) 

where at is the creation operator for a pion with momentum k n . The calculation of p Pd (k^) 
from a ttNN model will be given in the next section. 



IV. CALCULATIONS OF NUCLEON AND PION MOMENTUM DISTRIBU- 
TIONS 

We first describe a nuclear model within which the procedure for calculating the nucleon 
and pion momentum distributions in the rest frame of the deuteron is described. We then 
explain how the these distributions can be used to calculate the momentum distributions in 
a fast moving deuteron, which are the input to our calculations of Eqs. fj42|) and (j4"5|) . 



A. vriViV Model 



We follow the ttNN studies [29l . l37j to consider a nuclear model defined by the following 
Hamiltonian 

H = H + V NN + [ H nNN + H\ NN } , (45) 

where Ho is the sum of free energy operators for N and n, Vnn is a nucleon-nucleon potential, 
and HnNN defines the virtual N — > ttN transition 

H wNN = K NN (i), (46) 

i=l,2 

where i denotes the i-th nucleon. Starting with the standard pseudo-vector coupling, the 
vertex interaction takes the following familiar form (omitting the spin-isospin indices) 

Knn{i) = J dkdpidpl S(pi + k-p() [ \pik)F nNN (p i ,k)(p i '\ ], (47) 
where \pi) and \k) are the plane wave states of the z-th nucleon and pion, respectively, and 



(48) 

Here F(A n NN, k) is a form factor that satisfies F(A w nn, k) = 1 at k — in with k = 
m^^Jl — m^/(4m 2 N ) being the pion momentum at the nucleon pole of the ttN amplitude, 
and its cutoff parameter A n ^N can be determined in the fit to the tcN scattering data. 

It follows from Eqs. ( T45|) -(l471) that the bound state \^d) in the deuteron rest frame (pd = 0) 
is defined by 

H\V d ) = E d \V d } : (49) 
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where the deuteron wave function is normalized to ($d\^d) — 1 and has two components 

\^d) = t^[\(Pnn) + a\(p nNN )}. (50) 

Here Z is a normalization factor and each component of the wave function is normalized to 1: 
(4>nn\4>nn) = 1 and ((P-kNN^-kNn) = 1- By using the orthogonality condition (0 aw | 4*kNn) — 
0, Eqs. <M) and (EJJ lead to 



a(07rAw|(-ffo + VwaOI^ttAw) + (^wNnIH-^nnI^Nn) = dE d , (51) 

((p NN \(H + V NN )\(p NN ) + a((j) NN \Hl NN \(j) nNN ) = E d . (52) 



From Eq. (15"Tj) . we have 



-Ed — ((frwNN | (-Ho + VzViv) I^ttJVJv) 



(53) 



It is a difficult three-body problem to solve Eqs. (j5ip and f )52|) exactly and find a model 
of V/vat to fit the iViV scattering data and the deuteron bound state properties. Here we are 
simply guided by the results from the previous ttNN studies It was found that in 

the low energy region, the ttNN component is much weaker than the NN component, and 
it is a good approximation to neglect the matrix element of (^ttnnIVnnI^ttNn) in Eq. f[5"3"j) . 
We then have from Eq. (|52p . 



(4>NN 


H-kNN \4>-kNn)(4>ttNN 


H-itNNWNn) 


Ed — (<PnNN 


Ho\4>ttNn) 



((fiirNN 


HttNnI^nn) 


Ed — (4>ttNN 


Ho\4> n NN) 



\9NN\{-tia + Vnn)\(Pnn) H p 7- iu i j, \ — ^d, (54 J 

— {<PnN 

and the coefficient a of the total wave function ( 15U1) is 

(55) 

If we further assume that the pion loop (pion is emitted and absorbed by the same nucleon) 
in the second term of Eq. fj54j) can be absorbed in the physical nucleon mass, Eq. ( I54jl is 
equivalent to the following Schrodinger equation 

[H + V NN + VZ P {E d ) } \<f) NN } = E d \(f) NN ), (56) 

with the one-pion-exchange potential defined by 

~7/J — \(PtvNN\-Ho\(P-itNN) 

By using Eq. (HST) for H^^n, Eq. ( 1531) leads to 

\a\ 2 = (<l> NN \\f%<*{E d ) +p e * c (E d )]\(l )NN ), (58) 



where 



etc/n \ \ v 7 f /.x 7 rA r iv)(07rA r Af| , / .x / rri x 

Pr(^d) = 2>W*)7p 7T ^-rr r^h nNN (j), (59) 

%z r3 



i 



loop/77 \ u\ f\ \4>irNN)(4>nNN\ , , .x / ftn x 

An- P (^) = 2^ KnnWJp 77 VRU \s2 h ^N{t). (60) 
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FIG. 3. One-pion-exchange interaction in the center of mass frame of NN. 
From Eqs. ( 15 7p and ( I59p . we then have the following interesting relation 

PT(E d ) = -^VZ P (E d ). (61) 



Note that the relation Eq. (IBTI) is identical to that used in Ref. [38| to calculate the so-called 
pion-excess, except that a non-relativistic form of V^ p is used in their calculations. 

To see the physical meaning of p l ° op (E d ) and p e * c (E d ), we define the pion number N v in 
the deuteron rest frame as 



with 



From Eq. (150]) . we then get 



N w = J dkp„(k), (62) 



p n (k) = (y d \aiaz\V d ), (63) 



N n = ^\a\ 2 . (64) 

By using Eqs. ( 15"5I) and (I64"|) . we can define 

N^N^ + N^, (65) 

where 

K = ^(<f>™\(%° p {E d )\<l> NN ) 

dkpUk), (66) 

Nr = ^(<l>NN\pT(E d MNN) 

dkp e * c (k). (67) 



In the DY and DIS calculations, the contributions from p° n (k) are included in the meson cloud 
contributions to the nucleon parton distributions. Only p^ c (k) is needed in our calculation 
of pion contribution to proton-deuteron DY process. This assumption is similar to that 
used in the calculation of pion-excess contribution [38| to DIS cross sections 39| . 
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To calculate p^ c (k), we use Eq. f lBTj) by first calculating the matrix element of one-pion- 
exchange potential ( 1571) in the the rest frame of the deuteron. From the kinematics shown 
in Fig. [31 we have 



(<f>NN\VZ P (E d )\<pNN) = J dk(<t> NN \V$?(k,E d )\<j> NN ) (68) 
where the contribution from the pion with momentum k is 
((pNN\V^ p (k,E d )\(j) NN ) = [ dp0* NN (p + k)F* NN (p,k)- 



E d - E N (p) - E N {p + k) - E n (k) 
x^m(-P- k, k)<p NN (p). 

By using Eqs. ( ISTl) and (j48l) and including spin-isospin indices, we readily get 

= ^\~^ d ^NN\VZ P ^E d )\<p NN )} 

/ \ 2 

2 /"j- 1 ( Unn\ 1 m N m N 



dp- 



Z J (2vr) 3 \ J 2E„(k) E N (p) E N (p + k) 



E d - E N (p) - E N (p + k) - E w (k) 



x [F(A nNN , k)} 2 <p J N d ^(p, m s{ ,m s , 2 )(T d M d \I( n , r 2 )\T d M d ), (69) 

where the overall factor 2 comes from summing up two possible pion-exchange diagrams. 
The isospin matrix element is 

(Td,M d \I(Ti, r 2 )\T d M d ) = Yl Yl (TdM i \T a (l)\[Tr a N 1 N2] T , Mi ) 

T'M' T 

x(KiV 1 iV 2 ] r , M Jr_ Q (2)|T d M d )(-l) Q , 



a=-l,0,+l T'M' T 



with 



\TM T ) = (^rn l m 2 \TM T )\m 1 ,m 2 ), (70) 



v 22 

mi ,m2 



T'Mlj, J 



t 
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= ^ (^rnxm^tmt) (tlm t m a \T' M?) \m u m 2 ) , (71) 

t m\,m2 

where (jij2^i^2|j^) is the Clebsch- Gordon coefficient. 

For deuteron T = Mt = and only T' = M' T = for ttNN should be kept, we then 
get {T d M d \I(Ti,T 2 )\T d M d ) = —1/3. The spin-orbital part of the deuteron wave function in 
Eq. (169]) can be expanded as 

<pNN Jd (P, m si,™>8 2 ) = Y(LSM L m s \J d Mj d )(~m Sl m S2 \Sm s )u L (p)Y LML (p), (72) 



LM L 
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where the radial wave functions are normalized as 



p 2 dp[u 2 (p) + u\{p)} = 1. 



(73) 



Because the ttNN component is much smaller, the normalization factor is Z ~ 1. We will 

— * 

use Z = 1 for calculating p^ c (k) of Eq. (|69|) . In the same approximation, we will use the 
deuteron radial wave functions uo(p) and v,2(p) generated from the available realistic NN 
potentials such as ANL-V18 

Neglecting the small contribution from ttNN component, the nucleon momentum distri- 
bution in the rest frame of the deuteron can be written as 



p N (p) = ($ d |4V|$d) 



~ W Y A 



(4>NN\bp)?\<j)NN)- 



By using Eq. (!72|) and setting Z ~ 1, we obtain 

1 



Pn(p) 



An 



u o(p) + U 2(P)]- 



(74) 



(75) 



B. Momentum distributions in a moving deuteron 

In the calculation of proton-deuteron DY cross sections, the momentum distributions 
p Pd (f>N) in Eq. ( 1421) and pp^k^) in Eq. ( 1431) are defined in a fast moving deuteron with a 
momentum p^. To calculate such momentum distributions, we first note that the parti- 
cle number in a system is independent of the frame. We thus have the following frame 
independent normalization condition 



'ppAp 



dppAp) 



(76) 



where N a is the number of the considered particle a = N or tt in the deuteron, and the 
deuteron momenta (set pd in the z-direction) in the moving frame and the rest frame are, 
respectively, 

Pd = (E d (p d ),0,0,p z d ), (77) 
p° d = (m d , 0,0,0). (78) 

The nucleon momenta in Eq. (1761) are related by the Lorentz transformation defined by the 
velocity /3 = P z /E d (P d ) of the moving frame. Explicitly, we have 

E d (p d ) 



Pz 



E N (p) 



m d 
E d (p d ) 



Pz~ 



P d 

E d (p a 



E N {p) - 



Pz 



m d l E d (p d ) 
P± = P±- 

The above equations lead to the following Lorentz invariant relation 

; -»' r -> 



E N (f) E N (p)' 



(79) 

(80) 
(81) 

(82) 
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By using Eqs. ( ITS"]) and (|82|) . we thus have 

PpAv) = ^§)P^M (83) 

With Eqs. fl?9l - fl8"Tl) . we can use Eq. fl83|) to get p Pd {j?) from the momentum distribution 
p p °(p) in the rest frame of the deuteron; p p °(p) can be calculated from the momentum 

distributions in the deuteron rest frame: p e * c {k) of Eq. fl69l) for pions and Pn(p) of Eq. f!75|) 
for nucleons. 

Here we mention that the relation fl83l) for a two-nucleon system can be explicitly derived 
from the definition fl39|) within the relativistic quantum mechanics developed by Dirac, as 



reviewed in Ref. 40 



V. CALCULATION PROCEDURES 

In this section, we develop procedures to apply the formula presented in previous sections 
to investigate the Fermi motion and pion-exchange effects on the ratio R p d/ PP = <J pd / (2cr pp ) 
between the pd and and pp DY cross sections. Our first task is to relate our momentum 
variables p p , px and q to the variables used in the analysis [E| of the available data. This 
will be given in Sec. IV Al The procedures for calculating DY cross sections are given for pp 
in Sec. IV Bl and for pd in Sec. IV CI 



A. Kinematical variables for DY cross sections 

It is common [5] to use the collinear approximation to define the parton momentum: 

Pq P = XlPp , (84) 

P qT = x 2 p T , (85) 

where p qp (p qT ) is the momentum of a parton in the projectile (target), and x\ and x 2 are 
scalar numbers. The momentum q of the virtual photon in the qq — > 7 — > p + p~ , as seen in 
Fig.d is 

1 = Pq P + Pq T = X lPp + X 2PT- (86) 

In the considered very high energy region E p > 100 GeV, the masses of projectile (p^ = m^) 
and target (p? = can be neglected and hence s = (p p +Pt) 2 ~ %P P -PT, Pr-q ~ xiP p -Pt, 
and p p ■ q ~ x 2 p p • Pt- We thus have the following relations 

x, ~ (87) 

x 2 ~ (88) 

It is most convenient to perform calculations in terms of X\ and x 2 in the center of mass 
system in which the projectile is in the z direction and the target in —z direction: 



Pv = (Jp 2 + m2,0,0,p) ~ (p,0,0,p), (89) 



Pt = (\ P 2 + m 2 T , 0, 0, -p) ~ (p, 0, 0, -p). (90) 
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With the choices (1591 and (!90|) . we have 

s = (Pp+Pr) 2 ~ V, 
M 2 = q 2 = (xiPp + x 2 Pt) 2 ~ 4xix 2 p 2 . 

The above two equations lead the simple relation 



s 

By using Eqs. (I87I) - (I90I) . we can define a useful variable xp 

2(Pt-Pp) ■ q 

Xf = Xi — X2 ~ 

s 

■ q 2y/sp p ■ q 



M 2 

X\X 2 ~ • (91) 



(92) 



In the notation of Ref. [5j], we write 



XF ~ 7^2' (93) 

where pj| = p p ■ q is clearly the longitudinal component of the intermediate photon momentum 
with respect to the projectile in the center of mass frame. Experimentally, s, M, xp, 
and do I [dMdxp] are measured. With the relation ( 191]) . we certainly can determine the 
corresponding xi, x 2 and da/ {dx\dx2)- We thus will only give the expression of do j (dxidx 2 ) 
in the following subsections. 

B. Calculation of pp DY cross sections da pp jdx\dx 2 

We now note that with the simplifications used in defining the variable X\ and x 2 , as 
described above, the flux factor associated with Eq. ( )29|) become 1. Substituting Eq. ( !T5|) 
into Eq. ( )29|) . the DY cross section for p(p p ) + N(pn) — > /i + + A* - + X p + Xn with N = p,n 
is then calculated from 

daPN{ f q f N) = E / dpMUlM)fM + f q MflM)\ 

x^fe>6(q>-(p q + Pq y). (94) 

In the chosen center of mass frame, defined by Eqs. fl89|) and fl90|) . let us consider q in 
the target nucleon moving with a momentum p^ = (p z N ,0,0,p z N ). In the precise collinear 
approximation, only z-component of the q momentum is defined by p z N . As defined by 
Eq. ( 1881) . we thus write f>q = (pgj_,p|) where 

P\ = x 2 p z N , (95) 

and Pq± can be arbitrary. The integration over the q momentum distribution in the target 
N can then be written as 

/ dnfUn) = J **>firM, ( 96 ) 
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with 



fhM = Pn 
Similarly, we can define for the projectile prot< 



;on 



d P JIM 



By using Eqs. fl9"6T) and f[9"8"j) . Eq. (|9"11) can be written as 
da pN (p p ,p N ) 

g 



(97) 



(98) 



dq 2 



J dx x dx 2 [P p {x x )f N {x 2 ) + f^x x )f%{x 2 )\ 
Ana 2 



x ■ 



]S(q 2 ~ ( Pq + Pq) 2 ) ■ 



(99) 



9q 2 q 

Integrating the above equation over dq 2 , we then obtain an expression of the cross section 
in terms of x x and x 2 that are defined by experimental kinematics 

A-tvcJ 2 



da pN (p p ,p N ) 



dx x dx 2 



E 



%Pq +Pq) 



(100) 



which is the same as Eq. used in the analysis of Ref. |5|] since (p q + Pq) 2 = q 2 = M 2 for 
the partonic process qq — > 7. Therefore we identify f^(x), defined by Eq. ( 1971) . with PDFs 
of the parton model [also for f^(x)]. To compare with the results of Ref. [5(, we use PDFs 
of CETEQ5m [35j in our calculations of Eq. fllOOp . 

Equation (11001) for the pp then obviously takes the following form 



da pp (p p ,p N , 
dx\dx 2 



Aira 2 
9M 2 



(/;(xi)/;(x 2 ) + f^ Xl )f;(x 2 )) 



(101) 



C. Calculation of pd DY cross sections of da pd /(dx\dx 2 ) 

We first consider the contributions from the nucleon momentum distribution P Pc1 {pn) t° 
Eq. ( I42p for the proton-deuteron DY cross sections . With the simplifications used in defining 
the variable x x and x 2 , as described in the subsection V.A, the flux factor associated with 
Eqs. ( 112]) become 1. We thus only need to consider 

da pd {pp,Pd) ^47r a 



dq 2 



47T ' Gi /*/*/* 
Yl J ^ J ^ J d P N PP^PN)f} p {P q )fl N {Pq)S {q 2 - {Pq+Pq) 2 ) ■ 

(102) 



The above expression is for the contribution from an anti-quark in the nucleon N of the 
deuteron and a quark in the projectile proton. Other contributions have the similar expres- 
sions, just with different quark indices. 
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FIG. 4. (color online) Nucleon momentum distribution p Pd (p z ) = J dp±p Pd (p z ,p±) in a deuteron 
moving with momentum p^ in the z-direction. The deuteron wave function of Ref. 30fl is used. 



By using the definitions of the parton distributions, Eq. ( 196]) for f 



N 



and Eq. 

for fp{xi), Eq. f 1 1 2 j) can be written in terms of momentum fraction variable X\ for q in the 
projectile proton p and of the nucleon N in the deuteron. We then obtain 



da pd (p p ,p d ) 
dq 2 



E 



Ana 2 „ 2 



W Pd (^v)/ p 9 (a;i)/^(x^)(5 (q 2 - (p q + Pg) 2 ) . 

(103) 



Similar to the pp case, the deuteron momentum is chosen to be in the ^-direction: p^ = 
(PdiPd± = 0)- Before we proceed further, it is necessary to relate the momentum fraction 
in Eq. (j!03p to X2 which is determined by the experimental variables M, s and Xf through 
the relations: X1X2 = M 2 /s and Xf = x\ — x% = pl/(y/s/2). Since our derivation is based 
on the impulse approximation that the parton is emitted from the nucleon in the deuteron, 



it is appropriate to assume that the momentum of the emitted parton is p| 
pl ve is the averaged nucleon momentum in the deuteron defined by 



(p, 



z \2 
ave/ 



I P P M(P z ) 2 dP 
fPpM d P 



X2Pl ve , where 



(104) 



Note that p Pd (p) in the above equation is calculated from the nucleon momentum distribution 
Pn(p) [Eq. (1751) ] in the deuteron rest frame by using the relation Eq. flHBl . In Fig. HJ we 
show the dependence of the calculated p Pd (p z ) = J dpj_p Pd (p z , pj_) on a deuteron momentum 
Pd- As expected, we find that that p^ ve ~ p d /2 at each deuteron momentum. 

Changing the integration variable by = pl/p z N = X2Pl ve /p z N , we can write Eq. fllQ3p 
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as 



da pd (p p ,p d ) 4na 



dq 1 



9q 2 I dxi I dx2 I d PNP Pd (PN) i ^rfp(xi)f q N (x 2 p z ave /p z N ) 



Pn 



x8(q 2 - (p q +p q f) 



(105) 



Integrating over q 2 on both sides of the above equation, we then obtain an expression of the 
cross section in terms of x\ and x 2 which are defined by experimental kinematics 



da pd {p p ,p d ) 4vra 5 



dx\d%2 9q 2 q p 

where the q contribution is isolated in 



Pd , N ^2) = I ap N p Pd {pN)^rf 9 N {^P z a Jp z N )- 

Pn 



(106) 



(107) 



Within the parton model, we should only keep the contribution from f N (x2P ave /p z N ) with 
x 2PIvc/Pn — 1- The above equation can then be written as 



*LnM = / dp^f N {x 2 plJp z N )p Pd {p z N l 



with 



P Pd (P 



dpN±P Pd (PN,PN±)- 



(108) 



(109) 



The derivation of Eq. (11081) can be extended to have q in deuteron and q in the projectile 
proton. We finally obtain 



da pd (p 



pi 



dx\d%2 



9q 



2 i 
q 



fpM E f LnM + E F l 

,n( x z) 

N=p,n N=p,n 



• (no) 



We use the charge symmetry to calculate PDFs for the neutron from that of proton: f d — f p , 

frt = fpi fn = fpi fn = fp- Furthermore p Pd {p z N ) is the same for neutron and proton. 
Including the charges for u and d quarks appropriately, Eq. ( 11101) can be written as 



d<J pd (p p ,p d ) 4vra 2 



dx\dx 2 



9q 2 



+ f Lp^) 

[^pM + f Lp^ 



(111) 



The formula for calculating the contribution from pion momentum distribution can be 
derived by the similar procedure. We obtain 



da pd (p p ,p d ) _ 4?ra 2 



dx\dx2 



9q 2 



9* 



9' 



;ii2) 
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where f% w {x) is PDFs for the pion taken from Ref. 28J, and the convolution function for the 
pion is 



with 



Here the average pion momentum is defined by 

/ Ppd (kir)(k*) 2 dk n 



z \2 
ave/ 



(113) 



(114) 



(115) 



The pion momentum distribution p Pd (k n ) in the above equations is calculated from p e * c {k) 
of Eq. (1551) by using the relation Eq. 



VI. NUMERICAL RESULTS 

As discussed in Ref. [rj, the ratio d/u in the proton can be extracted from the data of 
the ratios between the pd and pp DY cross sections: 

R d = da pd (p,p d ) I / 2 <fr"(p,j»p) \ ^ 
pd/pp dx\dx 2 I \ dx\dx 2 J ' 

where X\ and X2 have been defined in Sec. IV Al We are interested in the effects of pion- 
exchange and nucleon Fermi motion on this ratio. The pp cross section da pp (p, p p ) / {dx\dx 2 ) 
can be calculated from Eq. (jlOip . The pd cross section da pd (p, pa) / (dxidx 2 ) is the sum of the 
nucleon contribution calculated from Eq. (11111) and the pion contribution from Eq. (11121) . 
To compare with the results of Ref. [BJ, the nucleon PDFs f p {x) of CETEQ5m [35[ is used 
in our calculations. The PDFs f%(x) for the pion is taken from Ref. [28| . 

From Eq. (11111) . it is clear that the nucleon Fermi motion effects are in F£ N (x2) defined by 
in Eq. (HDBI>. If we set p z a Jp z N -> 1 in Eq. §M§, F q pd N (x 2 ) f N {x 2 ) since J dp N p Pd (p N ) = 1 
as defined by the normalization of states. The calculation of Eq. (11111) with Fp d N (x 2 ) — > 
fli{x 2 ) is then identical to that based on Eq. of Ref. [1]. The differences between this 
calculation and that from using Eqs. (II lip and (1108P will indicate the importance of nucleon 
Fermi motion effect on pd DY cross sections. 

To calculate the pion contribution with Eq. (11121) . we need to first evaluate F p q dn (x 2 ) 

defined by Eq. (I113p . The pion momentum p Pd (k n ) in Eq. (11131) is calculated from using 
the relation Eq. (1831) and p e * c {k) defined by Eq. (169|) . We see from Eq. (l69l) that the pion 
momentum distribution p e * c {k) depends on the tiNN form factor F(A n NN, k) [Eq. (1481)]. 



Following the previous tiNN studies [29|, [37[ , this form factor must be consistent with it N 
scattering data. In this work, we apply the 7iN model formulated in Ref. jUJ to determine the 
F(A 7T nn, k) by fitting the nN partial wave amplitudes [42| up to invariant mass W = 1.3 



GeV. The nN scattering within this model has been given in Ref. [41| and will not be 



repeated here. Our fits are shown in Fig. [5j The resulting parameters are not relevant to 
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this work and are therefore not presented. For our calculation, we only need the resulting 
tcNN form factor. 

We see in Fig. that the resulting irNN form factor can be fitted by the following 
modified dipole form 



F(A nNN , k) = ( *\ NN f 2 ) [1 + a(l + k 2 /K 2 )] exp [-6(1 + k 2 / k 2 )] , (117) 

. A wNN + fc • 



where re = m^l - K/(4m^)], A nNN = 685.7 MeV, a = 1.67 x 1(T 3 , b = 2.79 x 1(T 4 . It 

/A 2 —k 2 \ 2 

is close to the usual dipole form F(A 7I -nn, k) = ( A 2 JVJV +fc2 J with A^^n = 810.6 MeV. 

The pion momentum distribution p^ c (k) calculated from Eq. ( 169]) with the ttNN form 
factor given in Eq. f)117p is the dashed curve in Fig. [71 Here we also show the nucleon 
momentum distribution Pn(p) (solid curve). Note that p c * c {p) changes s ign at p ~ 200 
MeV. This sign change is also seen in the calculation of pion-excess in Ref . 38( , except that 
their magnitudes are much larger because they use a much larger ttNN cutoff A ~ 1400 
MeV for a dipole form of a non-relativistic iViV potential. 

With the input specified above, we can calculate ratio R p d/ PP defined by Eq. ( 11 16j> . We 
compare three results: (1) No nucleon Fermi motion (FM) and no pion-exchange (7r-exc) 
from using Eq. (11111) with F£ N (x2) — > f^i^)] (2) With FM and no 7r-exc from using 
Eq. fillip ; (3) With FM and with 7r-exc from adding the results from using Eq. (II lip and 
Eq. dH2D- 
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FIG. 6. (color online) Filled squares are the ttNN form factor obtained from fitting the ttN partial 
wave amplitudes [42]. Solid and dashed curves are two different fits to the form factor using the 

parametrization of Eqs. (|117p and dipole form F(A n NN,k) = I \$ NN ,1.2 ) with A^nn = 810.6 
MeV 

In Fig. El the calculated R p d/ PP at 800 GeV are compared with the data of Ref. |5|]. Our 
results with no Fermi motion and no pion-exchange (dot-dashed curve) are similar to that 
presented in Ref. jEf. The differences between the dash-dotted and dashed curves are due 
to the Fermi motion of nucleon inside the deuteron. The solid curve also include the pion- 
exchange effects. All three results are close to the data. Clearly, the nucleon Fermi motion 
and pion-exchange effects are small in the region covered by this experiment. Our results 
shown in Fig. suggest that the simple formula Eq. ([2]) is valid to extract the d/u ratio in 
the proton in the small X2 < 0.3 region. 

To facilitate the analysis of the forthcoming data from Fermilab, we present our prediction 
at 120 GeV in Fig. |9j We see that the Fermi motion and pion-exchange effects are small in 
the £2 < 0.4. However these two effects are significant at larger x<i- We have observed that 
the rapidly raising effect due to pion-exchange is due to the fact that the parton distribution 
in the pion is much larger than that for the nucleon at large x, as seen in Fig. (TDJ Clearly, it 
is necessary to include the Fermi motion and pion-exchange effects to extract the ratio d/u 
in the proton from the data of R p d/ PP in the large X2 region. 

VII. SUMMARY 

For investigating the pion-exchange and nucleon Fermi-motion effects on the DY process 
in proton-deuteronQoc?) reactions, we have derived convolution formula starting with a nu- 
clear model within which the deuteron has NN and nNN components. The nucleon Fermi 
motion is included by the convolution of PDFs of the nucleon over the nucleon momentum 
distribution calculated from the NN component. The contribution from the ttNN com- 
ponent is expressed in terms of a convolution of PDFs of the pion over a pion momentum 
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FIG. 7. (color online) The momentum distribution A-Kp 2 p(p) of the pion (it) and the nucleon (TV) 
in the deuteron. Note that 47rp 2 p_/v(p) is multiplied by a factor 1/500. 

distribution that depends sensitively on the 7rNN form factor. With a ttNN form factor 
determined by fitting the ttN scattering data up to invariant mass W = 1.3 GeV, we find 
that the pion-exchange and nucleon Fermi-motion effects can change significantly the ratios 
between the proton-deuteron and proton-proton DY cross sections R p d/ pp = a pd /{2a pp ) in 
the region where the partons emitted from the target deuteron are in the Bjorken x<i > 0.4 
region. The calculated ratios R p d/ PP at 800 GeV agree with the available data. For analyzing 
the forthcoming data from Fermilab, we also have made predictions at 120 GeV. 
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